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Abstract. The Fermi surface topology of the organic superconductor A- 
(BETS)2GaCl4 has been determined using the Shubnikov-de Haas and magnetic break- 
down effects and angle-dependent magnetoresistance oscillations. The former exper- 
iments were carried out in pulsed fields of up to 60 T, whereas the latter employed 
quasistatic fields of up to 30 T. All of these data show that the Fermi-surface topol- 
ogy of A-(BETS)2GaCi4is very similar to that of the most heavily-studied organic 
superconductor, k-(BEDT-TTF)2Cu(NCS)2, except in one important respect; the in- 
terplane transfer integral in A-(BETS)2GaCl4 is a factor ~ 5 larger than that in k- 
(BEDT-TTF) 2 Cu(NCS) 2 . The increased three-dimensionality of A-(BETS) 2 GaCl 4 is 
manifested in radiofrequency penetration-depth measurements, which show a clear di- 
mensional crossover in the behaviour of H C 2(T). The radiofrequency measurements 
have also been used to extract the Labusch parameter determining the fluxoid inter- 
actions as a function of temperature, and to map the flux-lattice melting curve. 
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1. Introduction 



There is considerable current debate over the nature of superconductivity in quasi-two- 
dimensional (Q2D) crystalline organic metals JT], ^ ||. The most heavily studied 
members of this family of materials are the K-phase BEDT-TTF salts (e.g. k-(BEDT- 
TTF) 2 Cu(NCS)2) [0]. Whilst nuclear magnetic resonance ||, penetration-depth f||], 
tunnelling || and other experiments appear to suggest that the superconductivity 
in these salts may be d-wave-like and mediated by spin- density- wave- like fluctuations, 
some doubts have been cast by recent controversial specific heat measurements, which 
may suggest that the order parameter does not possess the required nodes ||. Several 
theories 0, [|], ||, || stress the importance of the details of the Fermi-surface topology 
in providing suitable prerequisites for superconductivity; if the Fermi-surface geometry 
and interactions are altered slightly, it appears that BCS-like s-wave superconductivity 
may be the dominant low-temperature groundstate |IJ . 

Clearly, it is of importance to study organic superconductors with slight variations 
in Fermi-surface topology so that the effect on the superconducting groundstate can 
be assessed. In this paper we therefore report magnetotransport and radiofrequency 
penetration-depth measurements of the superconductor A-(BETS)2GaCl4. Shubnikov- 
de Haas oscillations, magnetic breakdown and angle-dependent magnetoresistance 
oscillations indicate that the effective masses of A-(BETS)2GaCl4 and much of the 
topology of its Fermi surface are similar to those of the most heavily-studied BEDT- 
TTF superconductor, k-(BEDT-TTF)2Cu(NCS)2. However, the magnetoresistance 
close to 9 = 90° implies that the interplane transfer integral in A-(BETS)2GaCl4 is 
approximately five times larger than that in k-(BEDT-TTF) 2 Cu(NCS) 2 , suggesting that 
A-(BETS)2GaCl4 is less two-dimensional. This increased dimensionality is manifested 
in the superconducting properties of A-(BETS) 2 GaCl 4 ; a clear two-dimensional to three 
dimensional crossover is seen in the temperature dependence of H C 2. We have also 
used the radiofrequency measurements to extract various parameters related to the 
interactions between the fluxoids and to reveal the melting of the vortex solid. 



2. Background information 



When the innermost four sulphur atoms of BEDT-TTF JTJ are replaced by selenium 
to produce BETS (where BETS stands for bis(ethylenedithio)tetraselenafulvalene), 
electrocrystallization gives rise to a range of charge-transfer salts with notably different 
properties compared to their BEDT-TTF counterparts [TO, O, O, 113] . Salts of the A- 



phase morphology are presently unique to the BETS series, and A-(BETS) 2 GaCl 4 (T c 
5 K) remains the only superconducting BETS charge-transfer salt found thus far [[II 

Crystals of the A-phase exist in the form of needles with the long axis of the needle 



corresponding to the shortest lattice vector c pT|, |12l. At first sight, the crystal structure 
looks quasi-one-dimensional, with the BETS molecules packing roughly parallel in 
the planes between the anions [Q. However, the BETS sites are not equivalent, 
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and the cation molecules in fact occur in dimers, surrounded roughly isotropically 
by four nearest-neighbour dimers with the same orientation. The crystallographic 
unit cell contains two dimers, and hence contributes two holes | 12]| . Although the 



details of the cation positioning and symmetry are rather different than in k-(BEDT- 
TTF)2Cu(NCS)2 |Q, the overall similarity of the dimer arrangements lead one to expect 
a Fermi surface for A-(BETS)2GaCl4 which is topologically similar to that in k-(BEDT- 
TTF) 2 Cu(NCS) 2 . Indeed, the calculated bandstructure of A-(BETS) 2 GaCl 4 predicts a 
Fermi surface consisting of a quasi-two-dimensional (Q2D) hole pocket (the a pocket) 
and a pair of warped quasi-one- dimensional (Q1D) sheets |L2], [151 . According to the 
calculations, the a pocket in this case is expected to occupy ~ 28 — 33 % of the Brillouin 
zone 



12, 15 



3. Magnetotransport studies. 

3.1. Experimental details. 

Single crystals of approximately 1 x 0.1 x 0.05 mm 3 were synthesized using electro- 
chemical techniques |T2j employing a l,l,-trichloroethane/l,l,2-trichloroethane/ethanol 
solvent system fll6| . For the purpose of performing four- wire resistance measurements, 
12 jum gold leads were attached to the samples using graphite paint. In the pulsed field 
experiments, the resistance was measured using a 10 fiA ac current with a frequency 
of 200 kHz 0, |17|. The voltage was measured using a high-speed lock-in amplifier. 



Temperatures as low as ~ 340 mK were achieved by immersing the sample in liquid 



•^He inside a plastic cryostat fll7fl . Capacitor-driven, ~ 40 millisecond-duration pulsed 
magnetic fields of up to 60 T were provided by the National High Magnetic Field Lab- 
oratory (NHMFL), Los Alamos. Angle-dependent magnetoresistance (AMRO) studies 
were made using a two-axis rotation insert in quasistatic magnetic fields of up to 
30 T provided by NHMFL, Tallahassee. In the AMRO experiments, an ac current of 
5 fiA (frequency 30 — 80 Hz) was used for the resistance measurements, and a stable 
base temperature of 1.4 K was obtained by pumping on 4 He liquid. In both pulsed and 
quasistatic measurements, the current through the sample was driven in the interplane 
b* direction; in such a configuration the measured resistance is accurately proportional 
to the interplane resistivity component p zz |l|]. 

3.2. Pulsed-field magnetotransport: Shubnikov-de Haas and magnetic breakdown 
oscillations. 

Figure |I| shows the magnetoresistance of a A-(BETS)2GaCl4 crystal obtained using 
pulsed magnetic fields; the temperature was 340 mK and the magnetic field was applied 
parallel to the b* direction (i.e. perpendicular to the Q2D planes of the crystal) (Tl[ 



After the superconducting to normal transition, the resistance risies until a series of 
low-frequency Shubnikov-de Haas oscillations emerges at about 33 T. These grow in 
amplitude, until at about 45 T, a higher frequency series of oscillations becomes visible. 
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Field (tesla) 

Figure 1. Resistance of an A-(BETS)2GaCl4 crystal as a function of magnetic field, 
applied parallel to the b* direction; the temperature is 340 mK. The inset shows a 
Fourier transform of the data after subtraction of the non-oscillatory beackground |l[] . 
The peak at 650 T is associated with the a pocket; that at 4030 T is due to the 
breakdown (/3) orbit which encompasses 100 % of the Brillouin zone. 



The inset shows a Fourier transform of the magnetoresistance data. The lower of the 
two frequencies F a (believed to originate from the a pocket) is 650 ±5 T. The higher 
frequency of 4030 ± 25 T, which occurs at fields above ~ 45 T, corresponds to an area 
in &;-space approximately equal to the Brillouin-zone cross-section; following common 
usage in other charge-transfer salts Jl|, we will refer to this as the (3 frequency Fp. 
Frequencies equivalent to the Brillouin-zone area are readily observed in other charge- 
transfer salts (typically in k and a-phase salts of the form (BEDT-TTF^X) as a result 
of magnetic breakdown, whereby electrons tunnel between the Q1D and Q2D sections 
of the Fermi surface [I]. 

The observation of magnetic breakdown is not unexpected in this material, given 
the small size of the gap between the Q2D and Q1D Fermi-surface sections predicted 
by the bandstructure calculations [0. However, the experimentally-observed value of 
F a pa 650 T is roughly a factor two smaller than that predicted by the calculations. 
Discrepancies between model and experiment of this size are not unknown in charge- 
transfer salts |TR 
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Salt 


m* a /m e 


Fa (T) 


m*p/m e 




T c (K) 


Source 


re- (BEDT-TTF) 2 Cu(NCS) 2 


3.5 


600 


6.5 


3920 


10.4 


[21] 


A-(BETS) 2 GaCl 4 


3.6 


650 


6.3 


4030 


5 


present work 



Table 1. Comparison of magnetic quantum oscillation frequencies and effective masses 
in A-(BETS) 2 GaCl 4 and k-(BEDT-TTF) 2 Cu(NCS) 2 . The table shows the effective 
masses m* and and frequencies F a and Fp corresponding to the a and P orbits of 
the Fermi surface. 



3.3. Effective mass determination. 

Although the Fermi surface of A-(BETS)2GaCl4 is Q2D, the relatively small amplitudes 
of the oscillations in the magnetoresistance (see Figure [TJ) , and the absence of harmonics 
in the Fourier transform (inset to Figure |l|) imply that the Lifshitz-Kosevich (LK) 
theory should provide an accurate description of the temperature dependence of the 
oscillations Q. According to this theory, the thermal damping factor has the form 



where \ = 14.69 TK _1 and m e is the free electron mass. On fitting the amplitudes of 
the oscillations as a function of temperature (14 different temperatures ranging from 
340 niK to 3.0 K were used) we obtain the effective masses m* = 3.6 ± 0.1 m e for F a 
and trip = 6.3 ± 1 m c for Fp. 

Table [1] compares the effective masses and Fermi-surface areas obtained in A- 
(BETS) 2 GaCi4 and the most heavily-studied re-phase BEDT-TTF superconductor, re- 
(BEDT-TTF) 2 Cu(NCS) 2 H (see Section 3.2 of Reference for similar data on other 
re-phase BEDT-TTF salts). Note that the Fermi-surface parameters of the two salts are 
remarkably similar. Moreover, bandstructure calculations in both salts predict m* ~ 
m e |1^, |21|, whereas the observed masses are a factor ~ 3.5 bigger than this, indicating 
that interactions which renormalise the quasiparticle masses ]T], ^T], ^] are of similar 
importance in both materials. In both cases, the bandstructure calculations [O, EHI] 



predict that m 2m*, in reasonable agreement with the experimental values, and 
suggesting that the renormalising interactions influence both Q1D and Q2D Fermi- 
surface sections in a similar manner [[l], j2l| . The only marked difference between A- 
(BETS) 2 GaCl 4 and k-(BEDT-TTF) 2 Cu(NCS) 2 is the relatively high value of the Dingle 
Temperature in the former material; the rate of growth of the oscillations with increasing 
field (see Figure 1) suggests a Dingle temperature of Td ~ 3.2 ±0.1 K for the a pocket. 
Typical values of T D in k- (BEDT-TTF) 2 Cu(NCS) 2 crystals |], |3[] (and indeed in re- 



phase BETS salts [24]) are often a factor ~ 5 smaller than this, indicating that the 
impurity scattering rate in A-(BETS) 2 GaCl4 is relatively high. The reason for this 
difference is not yet clear. 
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Figure 2. (a) Magnetoresistance of an A-(BETS)2GaCl4 crystal as a function of 6, 
the angle between the applied magnetic field and b* . Data are shown for a number of 
different <fi angles (listed at right side of the Figure), where <f> is azimuthal angle between 
the plane of rotation and the b*c plane. The lowest three traces were recorded at 30 T; 
the rest of the data were acquired at 27 T. AMROs are observed as gentle oscillations 
of the resistance, periodic in tan0; the small peak at 90° is due to the presence of 
a small number of closed quasiparticle orbits on the warped Fermi-surface sections, 
(b) Illustration of the method of locating the AMRO resistance maxima (0 = 7° data 
from (a)). The slowly- varying background magnetoresistance has been fitted to a 
fourth-order polynomial in 6 and subtracted from the experimental data, leaving the 
oscillatory component. Peaks periodic in tanf? are plainly visible. 



3.4- Angle- dependent magnetoresistance oscillations (AMROs). 

Figure |](a) shows the magnetoresistance of a A-(BETS)2GaCl4 crystal as a function of 
8, the angle between the applied magnetic field and b*. Data are shown for a number of 
different <ft angles, where <fr is azimuthal angle between the plane of rotation and the b*c 
plane. Distinct AMROs are observed; their dependence suggests that they are caused 
by a Q2D Fermi-surface section |l| |19[ . In such a case, maxima in the magnetoresistance 
occur at angles 6i defined by [p], [L9|] 

6 , ib[|tan0 i = 7r(i±-)+A(0), (2) 

where i is an integer, k» is the maximum Fermi wave-vector projection on the plane of 
rotation of the field and b' is the effective interplane spacing (see Fig. 0(b)). On plotting 
the positions of the maxima 9i versus i, taking account of the correct sign of the f 



term (18|, 19], we obtain straight lines at all azimuthal angles in accordance with these 
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expectations. On choosing b' to be interlayer spacing (18.4 A) obtained from X-ray 
diffraction studies 0, [HJ, we obtain the locus for k\\ versus <p shown in Figure 0(a). 

A locus in the shape of a figure of eight is the usual result for a pocket of elliptical 
cross-section [[1], [18], For a pocket of ideal elliptical geometry, the locus of k\\ is given 
by 

fc|| = [^COS 2 (0 - £inc) + ^ S in 2 (0 - Cine)] 1 , (3) 

where £- mc is the inclination of the major axis of the ellipse with respect to b*c plane. 
The parameters k x = 4.86 ± 0.08 nm -1 , k y = 1.63 ± 0.01 nm" 1 and £; nc = 19 ± 5° yield 
the best fit (solid curves in Figure §(a)). 




Figure 3. (a) Locus of fc|| versus azimuthal angle (f> derived from fits of Equation g 
to the AMRO data. Data are points and the curve is a fit to Equation || with the 
parameters k x = 4.86 ± 0.08 nm" 1 , k y = 1.63 ± 0.01 nm" 1 , 6nc = 19 ± 5° and 
b' = 18.4 A. (b) The experimental cross-sectional shape and orientation of the Q2D 
Fermi-surface pocket, shown alongside (and to the same scale as) the Brillouin zone, 
reciprocal lattice vectors a* , c* and calculated Fermi surface of Reference jl5| . The 
pocket is described by the parameters k x = 2.43±0.04 nm -1 , k y — 0.815±0.005 nm" 1 , 
£inc = 19 ± 5° and b' = 2 x 18.4 A w 36.8 A(see text and Equation |). 



The values of k x and k y deduced from the AMRO data indicate an ellipse area 
corresponding to a Shubnikov-de Haas frequency of 2608 ± 60 T. This is almost exactly 
four times larger than the value of F a observed in A-(BETS) 2 GaCl 4 (see Table |l|). It is 
inconceivable that this value of ~ 2608 T could be the actual area of the a pocket, as 
it would then occupy ~ 66 % of the Brillouin zone. Quantum oscillation measurements 
remain the definitive method for obtaining Fermi-surface cross-section areas |], ^Ql • On 
the other hand, the AMRO measurements shown in Figure behave exactly as one 
would expect for a Q2D Fermi-surface pocket fl~9]| , with no evidence for any significant 
misalignment of the sample; i.e. the fits to Equation ^| are straight lines and the peak 
feature at 9 ~ 90° occurs at 90° for all azimuthal angles. 

A possible explanation is that the true interlayer spacing (as perceived by the 
quasiparticles) is double the unit-cell height, or that there is a modulation of the lattice 
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in the crystallographic b direction; this would result in an effective interlayer spacing of 
V = 2 x 18.4 A « 36.8 A. Such a modulation of the lattice could occur in the event of 
a charge-density-wave (CDW) or spin- density- wave (SDW) instability; however, there 
is as yet no other evidence for the presence of such a groundstate (c.f. numerous 
other charge-transfer salts in which CDWs or SDWs cause extensive modification of the 
quantum-oscillation spectrum ffl, If such a modulation exists, it is too weak to be 
picked up by a careful X-ray study at 115 K (Bruker-AXS SMART6000 CCD, complete 
sphere of data, sixty-second frames, 0.3 degree scans) [26|. However, this does not rule 



out the possibility that doubling occurs at a lower temperature p6| . 

Whatever the mechanism, a doubling of b' (see Equation ^) would result in the 
Fermi-surface parameters k x = 2.43 ± 0.04 nm _1 , k y = 0.815 ± 0.005 nm" 1 and 
£i nc = 19 ± 5°, yielding an ellipse area corresponding to a Shubnikov-de Haas frequency 
of 652 ± 15 T, in very good agreement with F a = 650 ± 5 T derived in Section 3.2. 
Figure |3|(b) shows an elliptical-cross-section Q2D pocket of this size and orientation 
alongside the most recent calculation of the Fermi surface |i~5"| , based on structural 



studies carried out at 17 K. The Q2D pocket measured experimentally is smaller and 
somewhat more elongated than that suggested by the calculations; it occupies ~ 16 % of 
the Brillouin zone, whereas the Q2D pocket of the calculation is 28 % of the Brillouin- 
zone area. However, such discrepancies between calculation and experiment are not 



without precedent in crystalline organic metals [18, 19 



3.5. Estimate of the interplane transfer integral 

We now turn to the small peak in the magnetoresistance component p zz observed at 
9 = 90° in Figure |[ Thus far, we have treated only the Fermi-surface cross-section in 
the a*c* plane (see Figure 0(b)). However, a small, but finite interlayer transfer integral 
will lead to a warping of the Fermi surface in the interlayer b* direction (see Section 2 of 
Reference |]IJ]). When the magnetic field is almost exactly in the plane of the warping, a 
few closed orbits become possible on the warped sections of the Fermi surface (e.g. on 
the "bellies" of the warped Q2D Fermi cylinders) 28), |3"0 . For certain orientations 



of an in-plane field, closed orbits will be possible on both the Q1D sheets and Q2D 
cylinders of a Fermi surface such as that shown in Figure §(b); at other orientations of 
an in-plane field, only the Q2D cylinders will be able to support closed orbits on their 
bellies. Such orbits are very effective at averaging the interplane velocity component v z , 



and hence lead to a peak in p zz p8| 



As the magnetic field is tilted away from the in-plane direction [9 = 90°), the closed 
orbits will cease to be possible when 6 = 90° ± A, where the angle A is given by 

v± 

A (in radians) ~ — , (4) 

where v± is the maximum interlayer quasiparticle velocity and vu is the intralayer 
component of the quasiparticle velocity in the plane of rotation of the magnetic 
field |28], ^9[ |30fl . If the quasiparticle dispersion E{k b ) in the interlayer direction is 
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assumed to follow a simple tight-binding model, E(k(,) = —2t±cos(k b b) |31|, where t± 
is the interlayer transfer integral, then 

v± = 2t ± b/h, (5) 

where we have used the relationship hv = V^Eik) ]3T| to obtain v±_ from E{k b ) |30[] . 
Equations [| and therefore show that there is a direct proportionality between 2 A, the 
angular width of the peak in p zz , and the interlayer transfer integral, t±. 



a 




Tilt angle 9 (degrees) 

Figure 4. Peaks in the interplane resistance R zz (proportional to p zz ) close to 9 = 90° 
in A-(BETS) 2 GaCl 4 (thick line; T = 1.4 K, B = 30 T; present work) and k-(BEDT- 
TTF) 2 Cu(NCS) 2 (fine line; T = 520 mK, B = 42 T; Reference @). In the case of 
k-(BEDT-TTF) 2 Cu(NCS)2, the rapid oscillations at the edges of the figure are angle- 
dependent magnetoresistance oscillations (AMROs). The fine lines superimposed on 
the A-(BETS)2GaCLt data show how the full width of the peak is defined. 



Figure |] shows how the angular width 2A of the peak in p zz is defined; the 
limits of the peak are defined by the intersections of extrapolations of the background 
magnetoresistance and the edges of the peak Similar data for k-(BEDT- 

TTF) 2 Cu(NCS) 2 from Reference |JD| are plotted for comparison; note how the peak at 
90° is much narrower. Given the similarity of their intralayer Fermi-surface properties 
(see Table 1), this comparison immediately suggests a much smaller t± in k-(BEDT- 
TTF) 2 Cu(NCS) 2 than in A-(BETS) 2 GaCl 4 . 

In order to obtain a quantitative estimate of t±, it is necessary to use reliable values 
of v\\ |28|, |30fl . Figure ||(b) suggests that when the in-plane magnetic field is close to 
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the a* direction, it is likely that both Q1D and Q2D Fermi-surface sections will be able 
to support closed orbits; conversely, when the magnetic field is well away from this 
orientation, only the Q2D cylinder will be able to support closed orbits. As most of the 
accurate information deduced from the experiments in the previous Sections concerns 
the Q2D Fermi-surface section, we concentrate on the range of <fr over which it alone 
will determine the peak at 9 = 90°. 

We assume an effective-mass-tensor approximation for the in-plane motion on the 
Q2D Fermi-surface section JTJ; 

h 2,2 7,2,2 

E = ?A + UbL (6) 

Here m\ and m 2 are effective masses for linear motion in the x and y directions. The 
cyclotron effective mass in such an approximation is m* = [mxm-i) 1 ^ [31]. Using the 



lengths of the axes of the elliptical cross-section of the Q2D Fermi-surface section derived 
from the AMROs (2.43 nm _1 and 0.815 nm" 1 ; see previous section- we have divided by 2 
to account for apparent doubling of the unit-cell height) and m* = 3.6 m c (Table [I]), we 
obtain ^ = C§^j and mim 2 = (3.6 m e ) 2 = 12.96 ml, yielding m\ = 10.73 m e , 
m 2 = 1.207 m c and Ep w 20.95 meV (c.f m\ = 10.59 m e , m% = 1.177 m e and 
E F « 18.4 meV in k-(BEDT-TTF) 2 Cu(NCS) 2 H). Hence, using ftv = V k £(k) g|, 
and constraining E = Ep, the velocities v\\(4>) may be derived for the Q2D Fermi-surface 
section. 

The widths 2A derived from data such as those in Figure ^ are plotted as a 
function of <fi in Figure |5|. The Figure also shows the prediction of Equation |] using 
v\\ derived from Equation § the only fit parameter is t± (see Equation |^). We have 
chosen to fit data for ranges of at which the Q2D pocket is expected to be the sole 
provider of closed orbits for in-plane magnetic fields; i.e. we avoided field orientations 
close to —a* (0 = 90°) and a* (0 = 270°) (dashed lines) at which the Q1D sheets 
might also be expected to provide closed orbits in an in-plane field (see Figure 0(b)). 
(Note that the experimental data show a strong peak close to = 90° and 270° 
(instead of the minimum predicted by the model) suggesting that the Q1D sections 
are indeed the dominant cause of the peak in p zz at these orientations.) The fit 
yields v± ~ 1200 ms _1 , so that t±_ ~ 0.21 meV ||33|| . A similar procedure has been 
carried out for k-(BEDT-TTF) 2 Cu(NCS) 2 , giving an interplane transfer integral of 
t± ~ 0.04 meV [|30[| . Therefore, the interplane transfer integral in A-(BETS) 2 GaCl4 is a 
factor ~ 5 bigger than that in k-(BEDT-TTF) 2 Cu(NCS) 2 . 



3.6. Summary. 

In summary, the magnetoresistance measurements indicate that the Fermi surface of A- 
(BETS) 2 GaCl4 bears a strong resemblance to that of k-(BEDT-TTF) 2 Cu(NCS) 2 within 
the highly-conducting Q2D planes. Moreover, the effective masses in the two salts are 
almost identical and the renormalising interactions are probably of a similar strength. 
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Azimuthal angle <j) (degrees) 

Figure 5. Full widths of peaks in p zz close to 6 = 90° for A-(BETS) 2 GaCl 4 (points) 
plotted as a function of <\>. The curve is given by Equations and |6| with 

t± = 0.21 meV. Note that we have chosen to fit data for ranges of 4> a t which the 
Q2D pocket is expected to be the sole provider of closed orbits for in-plane magnetic 
fields; i.e. we avoided field orientations close to —a* (<f> — 90°) and a* (<p — 270°) 
(dashed lines) at which the Q1D sheets might also be expected to provide closed orbits 
in an in-plane field (see Figure ||(b)) 



The width of the peak in the magnetoresistance close to 6 = 90° suggests that 
the interplane transfer integral in A-(BETS)2GaCl4 is approximately five times bigger 
than that in k-(BEDT-TTF) 2 Cu(NCS) 2 . This implies that A-(BETS) 2 GaCl 4 is a less 
two-dimensional material. 

4. Penetration-depth measurements. 

4-1- Experimental details. 

In the current experiments, the penetration depth was inferred by placing the 
superconducting sample in a small coil which is the inductive element of a resonant tank 
circuit |34], The exclusion of flux from the sample, and hence the coil, decreases the 
inductance of the circuit; consequently the resonant angular frequency, u = 1/V LC, 



will increase. The well-known properties of inductors lead to = where AA 4 



is the change in flux area, Aq is the area of the measurement coil, AL is the change in 
inductance and L is the total inductance. For small changes in inductance, = 2-^, 
where Af is the change in resonant frequency, fo is the initial resonant frequency, and L 
is the initial inductance. Through simple geometrical relations it can be shown that |37j] 

AA^ _ 2r s A\ - AX 2 

A c R 2 1 ] 
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where AA is the change in penetration depth, R is the effective radius of the coil, r s is 
the effective sample radius. Simple estimates for the sample sizes used in the current 
experiments show that the second-order term is negligible P7j, so that 

r s fo 

The A-(BETS)2GaCl4 samples used in the current study have a needle-like 
geometry [|T^]; a typical example had approximate dimensions 2 x 0.170 x 0.084 mm 3 . 
In order to maximize the filling factor and cross sectional area, the coil was made 
rectangular, with an effective area of 1.34 mm 2 ; for the sample mentioned above, the 
effective sample radius is half the shortest dimension (0.042 mm) (note that the long axis 
of the crystal is perpendicular to the coil axis). Calibration was achieved by placing a 
spherical superconducting sphere of known size in the coil ||. The sample was orientated 
in the coil so that the oscillating magnetic field was parallel to the crystallographic 
b* direction, i.e. perpendicular to the highly-conducting planes [11]. The tank-circuit 
capacitance was provided by a 30 pF mica capacitor, and the circuit was driven at 
/ ~ 25 MHz by a tunnel-diode oscillator |34, |38|]. The sample and coil were placed in 
a 3 He cryostat or dilution refrigerator. Quasistatic magnetic fields of up to 30 T were 
applied parallel to b*. 



80x10 
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Figure 6. Change in penetration depth versus temperature at zero field in A- 
(BETS^GaCU. Note that the cryostat used could only provide controlled, slow sweeps 
of temperature down to 4.2 K (data shown as continuous curve). Below this, data are 
recorded at fixed, stable temperatures (points). The error bars on the points give 
typical uncertainties, valid across the whole temperature range shown. 



Figure |6| shows AA of a A-(BETS) 2 GaCl 4 crystal (deduced from the frequency 
of the tank circuit using Equation H) as a function of temperature. The sharp rise in 
penetration between 4.5 K and 5.1 K gives a very clear indication of the superconducting 
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to normal transition; above this temperature AA ~ 80 fim (i.e. comparable to the 
shortest dimension of the sample), indicating that the radiofrequency fields penetrate 
the whole crystal once it is in the normal state. This is in agreement with estimates 
of the low-temperature normal-state conductivity for A-(BETS)2GaCl4 [0, which lead 
one to expect an in-plane skin depth 5 ~ 100 fim [B£ 



12. 



Measurements at low magnetic fields: determination of pinning parameters 



Wu and Sridhar [[|l[] have treated the repulsively interacting flux lines in a type-II 
superconductor as periodic, damped harmonic-oscillator potentials modulated by a rf 
field. The physical justification of their model is that the repulsive interaction of 
the fluxoids causes the Abrikosov lattice f42| , [43(1 to resist higher flux densities in a 
manner analogous to the way in which a two dimensional network of springs resists 
compression. In such a model, the Labusch pinning potential parameter a corresponds 
to the restoring force on fluxoids displaced slightly by the current density J induced 
by the radiofrequency field [4~I|: r/^- + ax = 4>qJ ■ Here x is the fluxoid displacement, 
f] is a damping parameter and <po is the flux quantum In other words, a small 

perturbation displaces the fluxoid from its equilibrium position against the restoring 
force provided by the repulsion from neighbouring fluxoids and the pinning potential. 
In the limit of small magnetic field (H -C B c2 ), the damping due to fluxoid viscous 



drag jy], |4J may be neglected, leading to a linear relationship between changes in 
the square of the penetration depth AA 2 and the magnetic induction B inside the 
sample fUR , 



AA 



-B(H). 



/i «(Tr— ■ (9) 

Figure [7| (inset) shows the field dependence of AA 2 at a temperature of T = 700 mK. The 
sample is well inside the mixed state for fields above /j,qH « 0.01 T, so that B = fioH; 
hence, the linear dependence predicted by Equation |9| fits the data well, yielding a 
Labusch parameter of a(700 mK) « 1.4 Nm~ 2 . 

Data similar to those in the inset of Figure |^ were acquired at a range of 
temperatures; the resulting values of a are plotted in Figure [7| (main figure) as a 
function of temperature. The solid line fit through the data in Figure [F| is the predicted 



temperature dependence from the two-fluid Gorter-Casmir |43|, |45| model, 

ji N 4" 



a oc 



(10) 



Figure |7| also highlights some of the differences between organic superconductors and 
the "high-T c " cuprates. As T — >• 0, the Labusch parameter in A-(BETS)2GaCl4 tends 
to 1.5 Nm -2 , almost four orders of magnitude smaller than that in YBa2CusO 



41 



Moreover, whereas the two-fluid model is able to describe the A-(BETS)2GaCl4 data to 
a fair degree of accuracy, the pinning-force parameter in YBa2Cus07 has been shown 
to follow the temperature dependence (1 — (T/T c 



\2\2 
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Figure 7. Inset: the determination of the Labush parameter a(T) by a straight line 
fit to low-field AA 2 vs. n^H data (T = 700 mK). Main figure: experimental values of 
a versus temperature (points); the curve is a fit to the two-fluid model expression. 



4-3. Measurements at high magnetic fields: deviation from Campbell penetration-depth 
behaviour and upper critical field. 

Figure |8] shows the measured change in penetration depth AA for A-(BETS)2GaCl4 at 
intermediate fields and a temperature of 700 mK. AA varies approximately as \/~H, 
as expected in the Campbell scenario for fluxoid motion |41], |46[. However, above a 
field which we label HqH*, AA deviates from the y/H dependence; this implies that the 
periodic harmonic well models are no longer applicable. At fields above [ioH*, 
AA follows the approximate field dependence H 2 , until the penetration depth saturates. 
The identification of the upper critical field in organic superconductors from 



conductivity data has been the subject of considerable debate [[48], [49(1 ; the transition is 
intrinsically broad, and phenomena such as a pronounced "hump" in the resistivity 
and negative magnetoresistance are observed close to H c2 [Q. (Note that similar 



complications also afflict the "High-T c " cuprates |50|, §TJ.) Recently, a concensus has 
emerged whereby most of the transition region between zero resistance and normal-state 
magnetoresistance is regarded as a property of the mixed phase (see References [E| 
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Field (T) 



Figure 8. The change in penetration depth signal vs. field for A-(BETS)2GaCl4at 
700 mK. The deviation from the Campbell penetration depth behavior is indicated as 
H* while the saturation of AA signal indicates H C 2 , determined from the straight line 
intercepts. 



and references therein); H c2 is then defined as the intersection of the extrapolations of 
the transition region and the normal-state magnetoresistance | 35| . 

In penetration depth measurements, the broadening is less severe since the 
measurement is not dependent on a macroscopic net current flow across the sample [41. 



49]. The pinned fluxoids probed by rf fields do not experience as large an electric field 



gradient and hence the dissipation associated with the normal core is reduced [41, H 



Nevertheless, the transition is still somewhat broadened, and so we follow the same spirit 
of the convention used in resistivity studies |35| and the GHz penetration depth studies 
of Reference |^9[ , defining H c2 as the intersection of extrapolations of the penetration- 
depth curves below and above the point at which the saturation occurs (see Figure |8]). 
We can be confident that the saturated behaviour is characteristic of the normal state, 
as it continues up to at least 30 T without further features. This strongly suggests that 
the whole of the sample is penetrated by the rf fields in the normal state, as suggested 



in Section 4.1 



The values of H C 2 and H* deduced from the pentration depth measurements are 
shown as a function of temperature in Figure §. Four different A-(BETS) 2 GaCl 4 samples, 
taken from four separate growth batches, were used in the study. There were negligible 
differences in their behaviour, and the data from the four samples in Figure |9| overlie 
each other, suggesting that the characteristic fields measured are intrinsic properties 
of A-(BETS)2GaCl4. We shall return to the temperature dependence of H c2 in a later 
section. 
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Temperature (K) 

Figure 9. Notional phase diagram of A-(BETS)2GaCl4 showing the upper critical field 
(triangles) and H* , marking the inflection points in the penetration depth (diamonds). 
The quasistatic field is applied parallel to b*, i.e. perpendicular to the quasi-two- 
dimensional planes of the sample. Points from four different samples are shown, often 
overlaying each other. The upper solid curve is H C 2 oc (T* — T) 1 / 2 ; the dashed curve 
is H C 2 oc (T C 2 — T). The lower solid curve is a fit of the two- fluid model expression for 
the flux-line lattice melting. 



4-4- Flux-lattice melting at H* 



We now turn to the change in behaviour which occurs at the field H* (see Figure 
We attribute the change at H* to flux-line lattice melting, as H* follows the (T c — T) 2 
dependence expected from the Gorter-Casimir two-fluid model . 

Additional support for this attribution comes from considering microscopic models 
of the melting process. Houghton et al. [52| have considered the elastic moduli of 



the flux-line lattice and proposed that melting occurs when the mean thermal flux- 
line displacement d(T) is a substantial fraction of the Abrikosov lattice parameter 



(20 o /v / 3fi) 1 



/2 



i.e. 



d(T) « c^i. Here cl is the Lindemann parameter [|53f| , a function 



used very generally in the description of solid-liquid transitions; typically cl 
The explicit expression for d(T) is 



d{T) 




t 



'4(V2-i; 



+ 1 



0.1-0.2. 



(i - 1) 1 /* (i - b) \(i-by/ 2 , 

where t = T/T c and b = B/B c2 . The parameter G t = (167rV(£;BT) 2 )/(0ij# c2 (O)) 
describes the importance of fluctuations in a given system; 7 is the anisotropy term 
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defined by 



7 



(12) 



with t; xy and £ 2 being the in-plane and interplane coherence lengths respectively, k 
is another Ginzburg-Landau parameter giving the ratio of penetration depth to the 

In this case we require k — j^. 

: 1.4 nm and 



1 yields £j 



coherence length for a particular orientation f43| , p2| . 
Using Ginzberg-Landau theory in the limit T —>■ 
A ^150 nm for A-(BETS) 2 GaCl 4 , resulting in k « 107. 

In this model, the melting field at each temperature can be interpreted as the point 
at which the product of cl and the lattice spacing is roughly equal to the average flux line 
displacement. Thus, an increasing field reduces the average inter- vortex spacing, thereby 
facilitating melting []5"2]] . Substituting fioH* « 1 T at T = 700 mK and k = 107 into the 
above equations yields a flux-line displacement of approximately 6.0 nm, roughly 12% 
of the Abrikosov lattice spacing. This implies that cl ~ 0.12, a value entirely typical of 
a solid-liquid transition [|3[ . 

In isotropic superconducting systems, the melting of the flux-line lattice occurs so 
close to H c2 as to be indistinguishable from it [Q. However, in materials such as A- 
(BETS) 2 GaCl 4 (and in k-(BEDT-TTF) 2 Cu(NCS) 2 ; see References @, 0), the large 
anisotropy of the superconducting properties permits the melting line to be observed 
over extended regions of the H — T phase diagram, well clear of H c2 . 



5. Discussion: comparison of A-(BETS) 2 GaCl4 with 
/t-(BEDT-TTF) 2 Cu(NCS) 2 ; dimensional cross-over. 

Having seen in Section |3] that the Fermi surfaces of A-(BETS) 2 GaCl4 and k-(BEDT- 
TTF) 2 Cu(NCS) 2 bear some striking similarities within the Q2D planes but have 
interplane transfer integrals differing by a factor ~ 5, it is interesting to compare their 
superconducting properties. 

work of Belin et al. E9 



The 



has shown that conventional resistivity 
measurements can yield unrepresentative values for the upper critical field of k- 
(BEDT-TTF) 2 Cu(NCS) 2 (see also |51], |6T|); the difficulties result from the dissipative 
mechanisms mentioned in Section O , which act to broaden the resistive transition [|1], 
Thermal conductivity, magnetisation and penetration depth measurements seem 
to be less susceptible to these problems and give a better reflection of the true H c2 |49[ . 
We have therefore compiled the H c2 (T) plot in Figure ITD] using available thermal 
conductivity [f|£j, magnetisation J56 , |57| and MHz [^5| and microwave (12 — 25 GHz) [fll| 
penetration measurements. There is some scatter amongst the data from different 
measurements, but all suggest that (dH c2 /dT) increases in magnitude as T increases, 
and in fact the power law H c2 oc (T c — T) 2 / 3 is quite successful in describing the data 



(Figure |10|). 

A material made up of weakly-coupled superconducting planes may transform 
from a three-dimensional system to what is in effect a series of two-dimensional 
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superconductors as the interlayer coherence length decreases with decreasing 
temperature The dimensional crossover occurs when the interplane coherence 

length £ 2 becomes shorter than the interplane spacing of the quasi-two-dimensional 
layers. Muon-spin rotation studies [[58] have shown that this transition occurs in k- 
(BEDT-TTF)2Cu(NCS)2 at magnetic fields ~ 7 mT, i.e. three orders of magnitude 
smaller than typical values of H c2 (Figure ID). This strongly suggests that the variation 
of H c2 shown in Figure |10| is typical of a quasi-two-dimensional superconductor consisting 
of weakly-coupled layers []56], |59], |62| . 




2 4 6 

Temperature (K) 
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Figure 10. Critical fields in k-(BEDT-TTF) 2 Cu(NCS) 2 , plotted on logarithmic 
(a) and linear (b) field scales. The data for comprise filled triangles (MHz 
penetration data obatined using the apparatus described in the current paper |55|| ) 
filled circles (microwave penetration studies; Reference [T^j; errors in /iqH C 2 values 
typically ±0.5 T), open circles (thermal conductivity data of Reference jf9|; errors not 
given) and shaded diamonds (magnetisation data; the lowest temperature point was 
determined by examining the attenuation of de Haas-van Alphen oscillations [p6| and 



the higher-temperature points are from the scaling studies in Reference 57 ). The solid 
curve is proportional to (T c — T) 2 / 3 , with T c = 9.1 K. The triangles are the irreversibility 
field from magnetisation [ |56| ; the filled squares and stars represent 2D melting from 
magnetometry and GHz studies [^9) (see also the NMR data of Reference |io| ) . The 
hollow squares are from muon-spin rotation pq] and denote the 3D-2D transition. 



Figure |9] shows that the upper critical field of A-(BETS) 2 GaCl 4 (conducting planes 
perpendicular to the applied magnetic field) has a linear region H c2 oc (T c — T) that 
spans from T c to approximately 1.9 K. Below 1.9 K, a definite change in the slope of 
the upper critical field occurs, and H c2 begins to follow the power law H c2 oc (T* — T)^, 
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with T* a fit parameter; powers ( in the range 0.5 — 0.7 provide an adequate fit to the 
data. 

The behaviour of H c2 in A-(BETS) 2 GaCi4 at temperatures below 1.9 K is therefore 
very similar to that of H c2 in k-(BEDT-TTF) 2 Cu(NCS) 2 over the whole temperature 
range shown in Figure [10], and is thus characteristic of a two-dimensional superconductor 
with weakly-coupled layers pDL . On the other hand, the linear variation of H c2 in 
A-(BETS) 2 GaCl 4 at higher temperatures follows the expectations of Ginzberg-Landau 
theory for three-dimensional superconductors |67|]- We therefore attribute the 
change in gradient at 1.9 K to dimensional cross-over from quasi-two-dimensional (low 
temperatures) to three dimensional (high temperatures). As we shall now show, this is 
entirely consistent with estimates of the interplane coherence length in A-(BETS)2GaCl4. 



The magnetic- field-orientation dependence of B c2 []37| in A-(BETS)2GaCl4 is 



qualitatively similar to the predictions of the Ginzberg-Landau anisotropic effective 



mass approximation 43, p3 



B c2 (9) = - B ^ 6 = 0) , (13) 
^cos 2 (#) +7~ 2 sin 2 (0) 

where 9 is the angle between b* and the applied magnetic field and 7 has been defined 
in Equation 0. Using Ginzberg-Landau theory , the in-plane coherence length may 
be estimated from the upper critical field when the magnetic field is parallel to b*: 

H c2 (9 = 0) = (14) 



Fits of the 9 dependence of H c2 at T = 1.75 K |54], £8| (the temperature at which 
the change in gradient in Figure || occurs) yield £ 2 ~ 1.85 nm, almost identical with the 
interplane spacing, and supporting our assertion that the change in gradient in Figure ^| 
at 1.9 K is associated with a dimensional crossover. 

Dimensional crossovers with the magnetic field applied perpendicular to the Q2D 
planes have been observed in artificial Q2D superconducting structures |)9L [70| and in 



organic superconductors such as k-(BEDT-TTF) 2 Cu(NCS) 2 [0, |59|; however, in the 
majority of these cases, the effect of the crossover is observed at magnetic fields less 
than H c2 . A-(BETS) 2 GaCl 4 is perhaps unique in providing the correct anisotropy for 
the crossover to be observed in the behaviour of H c2 (T). 



In Section |3.5| we demonstrated that the interplane transfer integral in A- 



(BETS) 2 GaCl 4 is a factor ~ 5 larger than that in k-(BEDT-TTF) 2 Cu(NCS) 2 . The 
greater "three dimensionality" of the bandstructure of A-(BETS) 2 GaCi4 compared to k- 
(BEDT-TTF) 2 Cu(NCS) 2 obviously manifests itself in the superconducting behaviour 
(compare Figures ^| and |10|); whereas A-(BETS) 2 GaCl4 exhibits 2D-3D dimensional 
crossover in its H c2 (T) behaviour, H c2 (T) in /t-(BEDT-TTF) 2 Cu(NCS) 2 is entirely 
characteristic of a Q2D superconductor. 
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6. Summary 

In summary, we have measured the Fermi surface topology of the organic superconductor 
A-(BETS) 2 GaCl 4 using Shubnikov-de Haas and angle-dependent magnetoresistance 
oscillations. The data show that the Fermi-surface topology of A-(BETS)2GaCl4 is very 
similar indeed to that of the most heavily-studied organic superconductor, k-(BEDT- 
TTF) 2 Cu(NCS) 2 , except in one important respect; the interplane transfer integral 
in A-(BETS) 2 GaCl 4 is a factor ~ 5 larger than that in k-(BEDT-TTF) 2 Cu(NCS) 2 . 
The increased three-dimensionality of A-(BETS) 2 GaCl4 is manifested in radiofrequency 
penetration-depth measurements, which show a clear dimensional crossover in the 
behaviour of H c2 . The radiofrequency measurements have also been used to extract the 
Labusch parameter determining the fluxoid interactions as a function of temperature, 
and to map the flux-lattice melting curve. 

We have observed a discrepancy between the angle-dependent magnetoresistance 
oscillation and Shubnikov-de Haas data which suggests that the true unit cell height at 
low temperatures is double that inferred from X-ray studies. At present, this has not 
been detected by other techniques. 

It is interesting to note that the anisotropies and Ginzberg-Landau parameters 
of the organic superconductors A-(BETS) 2 GaCl4 (this work) and k-(BEDT- 
TTF) 2 Cu(NCS) 2 @, |5|, |U span the typical values found in "High T c " cuprates 
such as YBCO and BISCCO 114311 . However, as the current work has shown, in 



contrast to the cuprates, the Fermi-surface topologies and complete phase diagrams of 
organic superconductors such as A-(BETS) 2 GaCl 4 and k-(BEDT-TTF) 2 Cu(NCS) 2 can 
be mapped out in detail using accessible laboratory fields. Moreover, details of the 
bandstructure in the organics can be related directly to the superconducting properties. 
The availability of a large number of organic superconductors of varying dimensionality 
and bandstructure fl], [23|] should potentially allow very stringent experimental tests of 
models of superconductivity in layered materials to be carried out. 
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8. Appendix: a note on the radiofrequency response. 

In the current paper we have extracted parameters which describe the superconducting 
state of A-(BETS) 2 GaCl 4 under the assumption that all of the apparent changes in the 
penetration of the radiofrequency field are due to the variation of A. It is therefore 
very important to assess whether this assumption is valid. Moreover, as little has been 
written in the literature about the radiofrequency techniques employed, it is useful to 
summarise the artefacts which can affect the experimental data. For future reference, 
we hope that it is also useful to provide estimates of some of the paramters used in the 



theory used to model experimental data |47 



Coffey and Clem [f|7]] have treated the behaviour of a superconductor in a rf field 
over a broad frequency range. Using their approach, the contributions to the penetration 
depth signal from surface-impedance and skin-depth effects can be evaluated as the 
temperature or applied magnetic field are varied. The model defines boundaries at 
which the surface-impedance and skin-depth effects become non- negligible; these are set 



by evaluating the flux creep factor j|7| 



The flux creep factor e is determined by u, the ratio of the fluxoid barrier height Uq 

2k B T- 



to the typical thermal energy k B T, v = The flux creep factor is then determined 



by: e — 1/ Io{v) 2 , where Iq is a zeroth-order modified Bessel function of the first kind |47J . 
e parameterises the degree to which thermal effects assist the motion of fluxoids. 

In the limit of large e {i.e. e ~ 1), thermal excitation causes the behaviour of 
fluxoids to approach that of completely unpinned fluxoids. In this case, the complex 
effective resistivity (p v (uj)) becomes a contributing factor to the measured change in 
penetration depth. p v {oj) |47]] is given by the expression 

e + (ujt) 2 + i(l — s)ujt 
1 + (ujt 

where P/ = B<f) /i] is the flux flow resistivity and ujt represents the product of the 
measurement frequency and the relaxation time of the normal-state quasiparticles. 

An order-of-magnitude estimate of Uq is given by considering the energy at 
which the harmonic-oscillator potentials of neighbouring fluxoids cross, yielding ffi 



Pv(u) = ; 2 Pf, (15) 



Uq ~ -^aL 3 . The charcteristic length L [^T| will be roughly equal to the Abrikosov 



lattice spacing, L = £, so that the pinning well barrier height Uq can be estimated using 

3/2 

(16) 



Equation [16| shows that U decreases as the field B increases, so that the surface- 
impedance and skin-depth effects will be most prominent at high magnetic fields. Within 
the superconducting state, the highest field {i.e. worst-case scenario) at which we 
make quantitative deductions about vortex behaviour is poH*, the field at which the 
penetration depth indicates a divergence from the Campbell regime. Substituting the 
value for 700 mK, we obtain U ~ 16 K. At this field and temperature, v ~ 4, leading 
to e ~ 0.005. 
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Equation [15] also shows that the value of ojt contributes to the complex resistivity. 
Taking a frequency oj/2ir w 25 MHz from the current experiments and r from 
measurements of the normal-state resistivity [12j (or the penetration depth: see 
Figure 0), we obtain ojt ~ 0.005. Therefore, the fact that both e and ojt are <C 1 



indicates that the surface-impedance and skin-depth effects have negligible impact [[17 
on the experiments on A-(BETS) 2 GaCl 4 reported in this work. 
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